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Abstract 



We study a system of N- bosons in the plane interacting with delta function potentials. After 
a coupling constant renormalization we show that the Hamiltonian defines a self-adjoint operator 
and obtain a lower bound for the energy. The same results hold if one includes a regular inter- 
particle potential. 

1 Introduction 

We consider N— bosons of mass m in R 2 interacting with delta function potentials of strength g. The 
Hamiltonian for the system is 



defined on Hn = the TV- fold symmetric tensor product of £ 2 (R 2 ) with itself. The problem is to 
make sense of this as a self adjoint operator. This is necessary in order that the global dynamics 
ipt — e~' lHt ij}o be well-defined. However the expression is quite a singular and one finds that the 
coupling constant g must be renormalized to have a chance of success. 

The problem is fairly well understood for N — 2. We give a treatment below which involves 
introducing a momentum cutoff, choosing a cutoff dependent coupling constant, and then showing 
that as the cutoff is removed the Hamiltonians have a self-adjoint limit in the sense of resolvent 
convergence. 

The N = 2 problem has also been considered by Albeverio, Gesztesy, Hoegh-Krohn, and Holden 
PP. They take a different approach which involves specifying boundary conditions when the points 
coincide. We show that our results are equivalent to theirs. 

Our main interest is in general N and the challenge is to incorporate the wisdom gained for the 
two particle case into the multi-particle setting. Our solution involves introducing fictitious particles 
known as angels which serve as markers for two-particle subsystems. This approach was developed by 
one of us in the papers [3], [Sj. The present paper is a rigorous version of this work. The main result is 
again a proof that the cutoff renormalized operators have a self-adjoint limit in the sense of resolvent 
convergence. We also obtain a lower bound for the Hamiltonian. 

This problem was previously considered by Dell' Antonio, Figari, Teta who also consider d = 3. 
Our results for d = 2 are very similar to theirs. However the proofs are rather different. They use 
a concept of T-convergence rather than resolvent convergence. Also their 'bare coupling constant' 
depends on N and momentum as well as the cutoff, whereas ours depends only on the cutoff. 




l<i<j<N 



(1) 
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The present work seems to have some advantages in simplicity and flexibility. As evidence of this 
we obtain the new result that essentially the same conclusions hold if we include a regular inter-particle 
potential in addition to the delta function. 



2 Two particles 
2.1 point interaction 

We begin with a discussion of the case N = 2. Taking mass m — 1 and passing to to center of mass 
coordinates we have the Hamiltonian 

H = -A-gS (2) 

on the space £ 2 (M 2 ). In momentum space 



m) ( P ) = p^(p) - J ^(iw ( 3 ) 

This operator does not map into £ 2 (R 2 ) and cannot determine a dynamics as such. 
Instead we consider approximate Hamiltonians 

(H A iP)(p) = P V(P) - MP) / PA(qM<l)dq (4) 

where p A is the characteristic function of \p\ < A. We define 

PfTp = f(f,ip) (5) 
(If ll/H = 1 this is the projection onto /.) Then we can write with H$ = p 2 



H A - H - J^y, P pa ( 6 ) 



This is a bounded perturbation of the self-adjoint operator Hq and and so is self-adjoint on D(Hq) 
(Kato's theorem pj). We define the resolvents 

R (E) = {H - E)- 1 R A (E) = (H A - E)- 1 (7) 

when they exist. If they exist as bounded operators one says that E is in the resolvent set of the 
operator. The resolvent set for Hq is C — [0, oo). Since the perturbation is rank one the resolvent 
Ra{E) can be explicitly calculated. For E e C — [0, oo) one finds that E is in the resolvent set for 
H\ — E if and only if 

(2tt) 2 ^V (pa,Ro(E) P a) (8) 

in which case 

R A (E) = ME) + L Va =i } = 7^ ? ) Pb. (e )pa (9) 
V (2tt) 2 3 a - (pa, R {E)p A ) ) 

Indeed if (JSJ) holds then an explicit calculation shows that the right side provides a bounded inverse 
for H A — E. On the other hand if (2n) 2 g^ 1 — (p A , R (E)p A ) then 

(H A - E)R (E)p A = (l- j^(p A ,R (E)p A )^j p A = (10) 

and Ro(E)p A ^ so E is an eigenvalue of H A — E and not in the resolvent set. 
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Now we introduce a new parameter fj, > and make the choice 

-l 



9A = <7a(/x) = (2tt) 2 ( f (p 2 + ^dp ] (11) 

Thus g>A goes to zero logarithmically as A — > oo. Then we can write 

R A (E) = R (E) + £ A (/A -E)- 1 P Ro(E)pA (12) 

where 

U(a,b)= f {p 2 + a)- l dp- [ {p 2 + b)- 1 dp (13) 

J|p|<A J\p\<A 



For a, 6 > we have 



A 2 A 2 
£ A (a, b) =7rlog(— + 1) - 7rlog(— + 1) 

= 7r l 0g (I + -l)- 7r l 0g (I + -l) 



(14) 



In the last step we have canceled the divergence in each term by adding and subtracting it log A 2 . Now 
it is a simple matter to take the limit A — > oo and get 

£(a,b)=wlog(b/a) (15) 

Theorem 1 

1. For E real and not in {— /j, 2 } U [0, oo) the strong limit R(E) = liniA^oo Ra(E) exists and is given 
by 

R(E) = R (E) + £(m 2 , -E)- 1 P Ue (16) 
where Qe € £ 2 (M 2 ) is defined by 

n E ( P ) = (p 2 - E)- 1 (i7) 

2. R(E) is invertible 

3. For E complex and not in {—/J 2 } U [0,oo) the limit R(E) = liniA^oo R\(E) exists. There is a 
self-adjoint operator H(p) such that R(E) = {H(p) — E)^ 1 . 



Proof. 

1. Under our hypotheses £(n 2 ,-E) = tt \og(-E/fi 2 ) ^ 0. Hence £a(-£, 7J 2 ) ^ for A suf- 
ficiently large and £(tj 2 , — i?) -1 = liniA^oo £a(a* 2 > — -E) -1 - We also have in £ 2 (R 2 ) the limit 
Qe = hniA^oo Ro{E)p\ The result follows. 

2. To show the null space of R{E) is {0} it is sufficient to find a dense set V C D(H\) such that 
for 77 £ V we have the existence of 77* = lim^^H^ — E)n. For then if R{E)\p = we have 

(77, = Jim A -oo((flA - E) v , Rk(E)ip) = { V *,R(E)iP) = (18) 

for all 77 £ T> and hence i[> = 



3 



For our domain T> we pick u e 5(M 2 ) so the the Fourier transform u is in C^(M. 2 — {0}). For u 
in this domain we have 

(H A u)(p) = p 2 u(p) - p A (p) J p A (q)u(q)dq (19) 

We argue that the second term converges to zero so that H A u — > Hqu. Since g A — > and 
||Pa|| = it suffices that / p A (q)u(q)dq = 0(A~ 1 ). 

To see this first replace p A (q) = pi(q/A) by p* A {q) = p*(q/A) where p* is smooth approximation 
to the characteristic function of the unit disc. The difference is 0(A~ n ) for any n, and so it 
suffices to show J p A (q)u(q)dq = C(A _1 ). 

Since u G C£°(R 2 — {0}) we have v(x) — \x\~ 2 u(x) in the same space and sou= —A q v for some 
v £ <S(R 2 ). Then after integrating by parts 

J p* A (q)u(q)dq = J (-A q p* A )(q)v(q)dq (20) 

This is 0(A~ 2 ) since |A g p^(g)| is (D(A~ 2 ) and v(q) is rapidly decreasing. 

3. This follows from the first two parts and a version of the Trotter-Kato theorem quoted in the 
Appendix. 



Remarks. 

1. The resolvent has a simple pole at E = —p? so H(p) has the eigenvalue — p 2 . The residue is 
the projection onto the eigenspace which we see is spanned by fi_„2(p) = (p 2 + p 2 )^ 1 . This is a 
bound state. 

2. Our approach to this problem follows a path well-known to physicists. The problem is usually 
cited as an example of dimensional transmutation in which a model without a length scale (the 
coupling constant g is dimensionless) upon renormalization gains a length scale (namely p^ 1 ) 

This phenomenon is expected to occur in gauge theories in four dimensions. 

3. Let us compare our result with the result of Albeverio et. al. [J. They consider —A on on 
£ 2 (K 2 \{0}) and then obtain various self-adjoint extensions by imposing boundary conditions at 
the origin. They obtain a family of self-adjoint operators indexed by a parameter a taking all 
real values. They also have an explicit formula for the resolvent (a Krein formula) which is just 
like our equation (|16|) except that instead of £,(p 2 , —E) = n\og(—E/p 2 ) they have the function 
(p. 99, equation (5.16)) 

4n 2 a — 2^(1) + it \og(—E/A) (21) 
Comparing we see that they agree exactly if the parameters are related by 

log/x = -2?ra + *(!) + log 2 (22) 



2.2 extension 

The previous results can be generalized to allow an additional potential besides the delta function. 
We consider 

H* = -A + v- gS (23) 

For simplicity we will assume v is a bounded function on M 2 . To define this we again start with 
approximate Hamiltonians in momentum space 

H* =H + v' -(27r)- 2 g A P PA (24) 
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where gx — <?a (/•*) is as before and v' = TvT 1 is a convolution operator (T = Fourier transform). 

I* 
'A 



Since = ||w|| = ||v||oo this is still a bounded perturbation and so is self-adjoint on D(Hq). 



Without the approximate delta function we have 

H 1 =H Q + v' (25) 

which is also self-adjoint on D(Ho) and satisfies and Hi > — ||f ||oo- 
Resolvents are denoted 

Rx(E) = (H, - E)- 1 R*{E) = (H* E)~ l (26) 

If E is complex and not in [— ||u||oo, oo) then E is in the resolvent set for H\. As before we find that 
that such E are also in the resolvent set for fff if and only if (27r) 2 g^ 1 ^ (p\, R\(E)piC) in which case 

R*(E) = f£i (-E) + ( - Pfl l(E)pA (27) 



(2^)^ L -{pA,Ri{E) PA ) 



Theorem 2 

J. Jbr reaZ E < —eg with 



the strong limit R#(E) = limA^oo R^{E) exists. 
2. R#(E) is invertible. 



e = maxflMloo + 1, /i 2 e l|l,|l ~ +1 ) (28) 



3. Rft(E) — limA^oo Rj[ (E) exists for all complex E not in [— e ,oo). There is a self-adjoint 
operator H#(p) satisfying H#(p) > — eo such that 

R*(E) = (H*(fi)-E)- 1 (29) 

Proof. In the denominator in (|27|l we insert 

R^E) = Ro(E) - R 1 (E)v'R (E) (30) 

and find 

R*(E) = Ri(E) + ( 1 ) P fll(B)pA (31) 

As A — > oo we have in £ 2 (R 2 ) 



lim i?i(S)p A = lim i? (^)pA - R 1 {E)v'R (E)p A 

A— »oq A— >oo 

=o B -Ei(£;yfi B (32) 

Thus we have the limit R#(E) = limA^oo R\(E) given by 

R#(E) = R 1 (E)+( 1 x W, a (33) 

\£(m , --E'J + ("Ei « "s) ™ (fijs, v'Ri{E)v'il E ) J 

provided the denominator does not vanish. However HHsH^ = tt | £7 1 1 and since E < — ||u||oo — 1 we 
have ||iii(J5)|| < 1 and hence 

m E ,v'n E )\ ^l^llvlU < 7t 

Thus we can avoid vanishing provided £(^ 2 , —E) > 7r(||v||oo + 1) or log(— 15 / p 2 ) > ||w||oo + 1- This is 
our condition — E > /i 2 e" 1 '" 00+1 . 

Thus part one is proved. The second and third parts follow as in the previous theorem. 
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3 Many particles 



3.1 bosons 

We now turn to the many particle problem. It is convenient to work with all possible values of N at 
the same time, even though the main interest is at fixed N. This means we are working on the Fock 
space H — ®^ =0 Hn- This has the usual creation and annihilation operators a*(f),a(f) defined for 
/ € C 2 (R 2 ). We also have a(p) = a(5(- — p)) defined on the domain V which is the dense subspace of 
H with only a finite number of entries and wave functions in the Schwartz space <S(R 2 ). For ip S T> 
the function p — ► a(p)tp is rapidly decreasing. (Note that a*(p) = a*(S(- — p)) is not an operator.) 

The Hamiltonian has the form H = Hq + Hj. The free Hamiltonian Ho is X^=iP? /2 on the Hn 
and is essentially self-adjoint on V n Hn- It can also be represented as a bilinear form on V x V as 1 

f P 2 
H = u)(p)a*(p)a(p)dp ui(p) = — (35) 

The interaction with interparticle potential —gS(x — y) is given in momentum space by the bilinear 
form on I? x T>: 

Hi = J a *(Pi) a *(P2) s (Pi +P2 -Pi - p' 2 )a(p 1 )a(p 2 ) dp 1 dp 2 dp' 1 dp 2 (36) 

However this is not an operator. 
To remedy this we introduce 

H A = H a + flj iA (37) 
For Hj t \ we add momentum cutoffs p\, take the coupling constant g\ = gh{p) as before, and define 

„ -9A f (P1-V2, ,P'\-P'l^ 

Hl > A J Pa{ ^^ )pa{ ^T- ] (38) 

a*(p' 1 )a*{p' 2 )5(pi +p 2 -p\ - p' 2 )a(pi)a(P2) dp 1 dp 2 dp' 1 dp' 2 

Changing variables to 

p = Pi+p 2 q= - (39) 

we find for the associated quadratic form 

( p p p p \ 

[ a ( 2 + l') a ( 2 ~ q '^> a( ~ 2 + q ^ 2 ~ ^ V dpdqdq> 
Applying the Schwarz inequality first in Fock space and then in the integral we find 

|(V,ff/,AVOI <2^i (/ J \H P - + q )a( P --qM 2 dpd q 



(40) 



5A 



2 



/ ||a(pi)a(p 2 )^|| 2 dpidp 2 (41) 



2(2tt) 

Here Ao = / a* (p)a(p)dp is the number operator. 

On the A-particle subspace Hat we have N — N and hence i?i,A is a bounded quadratic form. 
This determines a bounded self-adjoint operator on each Hn and hence Ha defines a self-adjoint 
operator on each Hn with domain D(H ) n Hat. Taking the direct sum we get a self-adjoint operator 
Ha on the full Fock space. 

lr rhis means (</>, Hotp) = J u)(p)(a(p)4>, a(p)ip)dp or as a quadratic form (?/>, Hqi/j) = J u>(p)\\a(p)ip\\ 2 dp 
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3.2 angels 

Next we introduce angels. We define 

H = C 2 (R 2 ) ® TL (42) 
which is Fock space with an angel. For / e £ 2 (E 2 ) we define %(/) : H — > 7Y and x*(/) : W — > 7t£ by 

X (/)(ft®V) 

These are creation and annihilation operators for angels, they are adjoint to each other, and they 
satisfy 

X(f)x*(h)=(j>h) 

X*(h) X (f)=h(f,-)®I 1 ' 

There is also the operator x(p) = xOK* ~ P)) defined say on the dense subspace V C 7i defined by 
V = S(R 2 ) <gs V. 

An equivalent representation is 

H = C 2 (M. 2 ,H) (45) 
Then I? is a subspace of S(M. 2 ,T>) and on this domain 

X (p)¥ = *(p) (46) 

Next we introduce: 
Definition 1 

BA = / PA ( Pl 2 ?2 ) x *( pi + P2) a (Pi) a (P2)dpidp 2 (47) 

Then Ba is an operator from Ti to 7i, and the key point is that it provides a square root for Hj A . 
Lemma 1 For A < oo 

1. B A defines a bounded operator on each subspace Hn ■ 

2. For ip *E T> we have in the representation \^5\j 

(B K 1>M = / M?) o(| + <?M~ - «ty <fc (48) 

5. On eac/i Hat: 

- 9 aB* a B a = H IA (49) 

Proof. The expression is naturally defined as a bilinear form. For ip G D and f £ 2) we have in the 
representation l|45|) 

(*,B A ^) = y PA(^y^)(*(pi +P2),a(pi)a(p 2 )V')^idp2 (50) 

Applying the Schwarz inequality twice we have 

1/2 / „ x 1/2 



\(%B A il>)\ < (J \ PA (i^±)\'\\q,( Pl + p 2 )\\ 2 d Pl dp 2 J (J \\a(p 1 )a(p 2 W\\ 2 dp 1 dp 2 j ^ ^ 

<||pa|| 2 ||*||||JVo^|| 
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Now specialize to if) € TLn and we see that B\ is a bounded bilinear form and hence a bounded 
operator. This establishes the first point. 
Next change variables in l|5U|) and obtain 

(*, B A iJj) = J PA (q) (*(p), a(| + g)a(| - g)V>) dp*? (52) 

which establishes (|48H . 

For (|49() it suffices to establish the identity as a quadratic form on T>. Inserting the representation 
of Ij48(l into — g\\\B A ip\\ 2 we obtain the representation (|40J) of (ip, Hj yA tp). This completes the proof. 

For later reference we consider the case A = oo with the operator 

B = rx, n \ i X*(pi + P2)a{pi)a(p2)dpidp 2 (53) 



V2(2ir) 

Lemma 2 B defines an (unbounded) operator on H.n D D(Hq) which satisfies for some constant C: 

||B^||<C||(flo+JV)^|| (54) 

For ip in this domain 

lim B A ijj = Btp (55) 

A^oo 



x / (w(pi) + l)(w(p2) + l)\\a(pi)a( P2 )^\\ dpidp: 



Proof. All the above representations still hold for ip e V, "F G t>. But now instead of (|51|l we have: 

< (/ (^(pi) + i)- i (^(P2)+i)- i ii*(pi+P2)irdpi^2) 7 

1/2 

(56) 

1/2 

< (j ( W (| + 9 + + g ) + iri*(p)||^ g j ik^o + Wll 

<C||*||||(H + iVo)V|| 

Here C = (f(u>(q) + l) -2 ^) 1 / 2 and in the last step we use the Schwarz inequality in q. This shows 
that B defines an operator on V D Hn satisfying the inequality l|54|l . Since T> n 7i n is a core for fZo 
on 7Ytv we can extend the domain to D(Hq) D 7ijv- 

For the second point we estimate K'F, (£? — B A )ip)\ as above. In the last integral over q we are now 
restricting to \q\ > A. Break this into two terms using 

{<! ■■ \q\ > A} C |g : |f +q\ > -} U jg : || - g| > || (57) 

With <5C*a = (/ 9>A (w(g) + l)" 2 d<?) 1/2 and V € 2? n H N this leads 

|(*,(B A -B)^)| < v / 2C A q|*||||(#o + A0V'll (58) 
This estimate extends to ip € D(Hq) n 7Yat. Then as A — > we have 5Ca — > and ||(Sa — B)ip\\ — > 0. 
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3.3 resolvents 

We return to A < oo and work out some consequences of the identity H49fl for resolvents. We define 

R {E) = {H - E)- 1 R A {E) = (H A -E)- 1 (59) 

These exist for ImE ^ and Rq(E) exists for E < 0. We want to find real E such that R A (E) exists 
as a means to isolate the spectrum of H A . 

To this end we also introduce the operators on Ti Ti 

H A (E) = ( Ho ~ a E *\ ) R A (E) = H^E)- 1 (60) 

Since B A is a bounded operator from Hn to TLm-2 we have that H A (E) preserves the subspace 
Hn © Hn-2- More precisely it is defined on (D(H ) n TLn) ®Hn-2 and is self-adjoint there. 

Lemma 3 For E < 0, R A (E) exists in B(TLn) iff R a (E) exists in B(H.n © Hn—2) in which case 

R.(E)-( Ra{E) -9aRa(E)B a \ 

M ' ~ \ -g A B A R A {E) g A +g 2 A B A R A (E)B* A J [bL) 

Proof. We omit the subscript A for the proof. First assume that R(E) exists. Then it is self-adjoint 
and has the form 

R(E) = ( « £ ) (62) 

for bounded a, P,5 and a, 8 self-adjoint. The statement that is the inverse says that a, (3* map into 
the domain of Hq and that 

(H - E)a + B*/3 =1 
(H -E)P* +B*5 =0 

Ba + g- x {3=0 ^ 
Bf3* +g~ 1 6 =1 

We ignore the second equation. The third equation says 

P = -gBa P* = -gaB* (64) 

Inserting the expression for p into the first equation and using —gB*B = Hi we get (H — E)a = I. 
Hence R(E) exists and equals a. Inserting the expression for P* into the last equation gives 5 = 
g + g 2 gBR(E)B* 

On the other hand if R(E) exists one can check directly that (|61|l provides a bounded inverse. This 
completes the proof. 

Now we give another version. 

Definition 2 For E < define a bounded operator on each TL^ by 

$ A (E) = gl l - B A Ro(E)Bl (65) 

Lemma 4 For E < 0, R A (E) exists in B(TLn ®H-n-2) iff ®a{E)~ 1 exists in B(Hn-2) in which case 

~ ( R () (E)+R (E)B* A ^ A (E)- 1 B A R (E) —R (E)B A $ A (E)~ 1 \ 

R ^ E >-{ -<S> A (E)-iB A R (E) * A (£)"i ) ^ 
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Proof. Again suppose that R{E) exists so we must solve the equations <|63[) again. This time we 
ignore the third equation. Then the second equation says that 

/3* = -R (E)B*5 (3 = -5BR {E) (67) 

Substituting (3* into the fourth equation gives (—BR (E)B*+g~ 1 )6 = I or $(E)6 = I. Hence ^(E)- 1 
cxists and equals 5. Substituting f3 into the first equation gives (Ho — E)a — B* $(£ , )~ 1 B Rq(E) = I 
whence a = R (E) + R (E)B*^(E)- 1 BR Q (E). 

On the other hand if $(£')~ 1 exists one can check directly that l|66|) provides a bounded inverse. 
This completes the proof. 

Comparing these results we have: 
Lemma 5 For E < 0, R A (E) exists in B(Hn) iff ^a(E)^ 1 exists in B(H,n-2) in which case 

R A (E) =R a {E) + R Q {E)BX^ A (E)- 1 B A R (E) 

i o (68) 

^aOE)- 1 =g A + g 2 A B A R A (E)B* A 



3.4 renormalization 

In view of the last result we can study the resolvent R A (E) on Hn by studying the operator $a(-E) 
on TLn-2- The advantage of this operator is that it can be more easily renormalized. 

First we Wick order moving creation operators to the left and annihilation operators to the right 
using [a(p), a*(p')] = S(p — p') and 

(H Q -E)- 1 a*(p)=a*{p)(H Q +Lu(p)~E)- 1 (69) 

The resulting identity is formal but a rigorous version can be had by regularizing a(p) — > a(5 K (- — p)) 
with approximate delta functions 5 K . We find 

$a(E) = $ ,a(£) + $iA E ) ( 7 °) 

where 



$0,a( E ) =9k ~ 2(27r) 2 / dpidp2 x *( pi + P 2 ^P 1 +P2)pa( P1 2 P2 f 



(71) 



2(2tt) 2 J mz * ' ' """" 2 ' Ho + ^+uj-2-E 
*lAE) =- 2(^)2 / d Pl dp 2 dpW2X*{pi +P2)x(p'i +P2)pa(^ T ^)pa(^ T ^) 
a*(p'i)a*(p' 2 )—— ■ — — — ; -a(p 1 )a(p 2 ) 

Hq + Ll>i + Ll>2 + Ll>[ + LU 2 ~ E 

4 

+6(pi ~ P'i)a* (p' 2 ) M — — —7 ^a(p 2 ) 

Here ct>i = u>(pi) = Pi/2, etc. These are bilinear forms on 2? x T>. By the methods of section 3.2 they 
determine bounded operators on each Hn for A < oo. But now we want to work uniformly in A and 
also include A = oo. 

To cancel the divergence in <& ,a{E) we change variables and write 

$ ,a(£) = (2^ 2 ( / (q 2 +MT 1 - / dpdq X *(p)x(p) n T 2/! , 2 f ) (72) 
\yi«|<A i| q |<A H +p 2 /4: + q 2 - E I 
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In the representation Ti — C 2 (R 2 ,H) this is 2 

MsvM-vr'lf tf+^>-'-/ * w/ ; + ,,_ E l »w (73) 



-(27r)- 2 a(M 2 ,^o+p 2 /4-^) *(p) 
As noted in (|14fl . £a has no divergence and we can define for A = oo: 

($o(£)*)(p) -(27r)- 2 e(^ 2 ,i/o+p 2 /4-^) *(p) 

^ x-l, /#o+P 2 /4-F\ T , v ( 74 ) 
=(4tt) Mogf ^ J*(p) 

Lemma 6 For F < — // , ^(F) is essentially self-adjoint on T>C\TIn and for \& in this domain we 
have 

lim $ a(£)* = $o(F)* (75) 

A^co ' 

Proof. For the essential self-adjointness it suffices to show that the domain contains a dense set of 
analytic vectors. (Nelson's theorem, 8 ). For analytic vectors we can take wavefunctions with compact 
support. 

The convergence is straightforward. One can use the inequality 

log(tf +P 2 /4 - E + p) - log(F + P 2 /4 - E)j *(p)|| 
< A- 2 ||(H +P 2 /4 - B)- 1 *^)!! < 0(A- 2 )||vI/(p)|| 
which follows using the spectral theorem. 

Next we work on ^/^(F). For A = oo it is defined without the p\ and denoted $j(J5). 
Lemma 7 For F < -1 and A < oo and * e X>: 

|(*,$/,a(F)*)| < 2(*,A 2 *) (77) 
Thus $j a(-E)) ^i(E) define bounded operators on TLjq and for $ £ Un Hat: 

lim $/,a(F)* = $/(F)* (78) 

A^oo 

Proof. We take ^/ ^(F) = $j 2 A (F) + $j 4 A (F) where the superscript indicates the number of creation 
and annihilation operators. For the first we have 

/ d Pidp2dp' 2 \ (a{p' 2 )^( Pl + P2 ),-— 1 , ^a(p 2 )t>( Pl +p' 2 ) 

2ir z J \ H + uii + oj 2 + ui' 2 — F 

/ d Pi d P2dp' 2 \\a{p' 2 )^(pi+p 2 )\\ \ ||a(Pa)^(Pi +^)H (79) 



/ dpad^||a(^)^H 1 \\a( P2 )*\ 

J UJ 2 + UJ 2 + 1 



<||7V 1/2 *|| 2 



2 In general if T = f \* (p)T(p)x(p)d-P defines an operator on H = £ 2 (R 2 ) ® W, then in the representation H 
C 2 (R 2 ,H) we have (T*)(p) = Tfoj)*(p). 
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Here in the last step we use the fact, noted in [2], that for h, h' £ £ 2 (R 2 ) and any c > 0: 

hip) 2—^- h'(q)dpdq\ < 7r 2 \\hh\\h% (80) 
p + q + C 

For the convergence we proceed differently. We use the estimate for e > 

Ipa(^^)pa(^ T ^) - 1| < 0{K- f -){^ + u 2 +J 2 + If (81) 

Then for $i,f 2 6Pwe have 

\(* u ($?l(E)-<S>?\E))* 2 )\ 
<0(A~ e ) [ dp 1 dp 2 dp' 2 \\a(p / 2 )^ 1 {p 1 +p 2 )\\- — 1 , , lU _J fl(pa)*a(pi+j£)ll 

J (LJl + UJ 2 + Ll> 2 + 1 j 1 



<0{/T<) / d Pl dp 2 dp' 2 - j—^Ha^)^^ +P2 )\\ ^ 2 + y 2 ||a(p 2 )^ 2 ( Pl +^)H 



(82) 



•(w 2 + i)*- e/2 K + i) 

<0(A- 6 )||JV '*i||||(iro + JVo)'*2||| 

where the last step follows by the Schwarz inequality. Specializing to T> n 7i at the estimate is uniform 
in H^ill = 1 and yields the convergence \\(&p A (E) — $^(£')) , J' 2 || — ► (In fact strong convergence 
holds since we have a uniform bound on the norms). 
For the second term we define 

f(p>,q> 7 p) = \\a(^+ q i)a(^- q i)*(p)\\ (83) 



and find 



(£?)*) < y dp 1 d P2 dp' 1 dp' 2 



\\a{p' 1 )a{p' 2 )^{p 1 +P2)\\— — —7- , ||a(pi)o(p 2 )^(p , 1 + p' 2 )|| 

Wl + CJ 2 + U\ + OJ 2 + 1 

If 1 

-8^2 y dpdqdp ' dq ' /(p , ^ , ^) g2 + (g , )2 + 1 /(P' g»p) 

<~ / dpij/||/(p',-,j>)||2||/(p J -y)||a 



(84) 



<^ll/ll2<^o*|| 2 

Again we have used JSOJ- This completes the bound, and the convergence follows by an estimate 
similar to (|8*2"|) 

To combine these we have : 
Lemma 8 

1. For E < —X, <$>(E) is essentially self-adjoint on T>C\Ti.N and for ^ in this domain 

lim $a(£)* = (85) 

A— >oo 
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2. Let E < — ejv where 

e N =max(l,/i 2 e 167rAr2 ) (86) 
Then for A sufficiently large or A = oo u>e /iai>e i/iai $a(£-) is strictly positive and for \P S 7ijv 

lim ^a^)" 1 * = ^(E)- 1 ^ (87) 

A^oo 

Proof. = $o(-E) + $r(-E0 is the sum of a essentially self adjoint operator and a bounded 

operator. The essential self-adjointness again follows by Kato's theorem. The convergence follows 
from our results (23, (JZBJ- 

For the second part under our assumptions £(/i 2 ,— E) — it log(— E//i 2 ) > 16ir 2 N 2 . Then since 
£,a(h 2 ,—E) converges to £(fJ. 2 , —E) we have for A sufficiently large (depending on E, fi) £,a(h 2 , —E) > 
12ir 2 N 2 . Since £a(o, b) is increasing in b we have for 'J' S 2? n 7i n ■ 

(*,$o,a(-B)*) > (2^)- 2 a(M 2 ,-£)ll*H 2 > 3iV 2 ||*|| 2 (88) 
Combining this with the bound |(^, $/ i a(£') , J , )| < 2Af 2 ||\f r || 2 we have for A sufficiently large or A = oo: 

(W,$a(-E)*) > ^ 2 ||*|| 2 (89) 

This gives the positivity and shows that $\(E) has a bounded inverse. Convergence on the core VCXHm 
for Q(E) and the uniform bound ||$a(-E') _1 || < A^ 2 imply the strong convergence for <I>a(£ , ) _1 . (See 
for example [S], p. 429) 
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3.5 resolvent convergence 

Now we can prove the main result (c.f. DelP Antonio, Figari,Teta |2] ) 
Theorem 3 

1. For real E < — ejv and ip £ Hn the limit R(E)ip = limA-*oo Ra(E)iP exists and is equal to 

R{E) = R {E) + R {E)B*<P{E)- 1 BR (E) (90) 

2. R(E) is invertible. 

3. For E complex and not in [— cn,oo) the limit R(E)iJj = liniA^oo Rj\{E)ip exists. There is a 
self-adjoint operator H(p) with > — epj so R(E) — (H(n) — E) . 

Proof. 

1. By lemma[8]if E < — ejy and A is sufficiently large then <J>a(-E) -1 exists as a bounded operator 
on Hn-2- By lemma|21it follows that all such real E are in the resolvent set of Ha on TCn and 

R A (E) = Rq(E) + R q (E)B* a ^a(E)- 1 BaRo(E) (91) 

We claim that BaRq{E) converges in norm to BRo(E). By the resolvent identity it suffices to 
prove this for any E < and we take E — —N and show Ba(Hq + N)^ 1 converges in norm to 
B(H + N)~ l . This follows by ffify. Taking adjoints we have that R (E)B A converges m norm 
to Ro(E)B* . We also know by lemma[H]that $a(-E) _1 converges strongly to $(£') _1 . Combining 
these results we have that R\(E) converges strongly to R(E) given by l|9U|l . 

2. As in the proof of theorem it suffices to find a dense domain of vectors ip £ Hn so that 
H\ijj converges. In fact we show Hi^ip — > which suffices. We have Hi^ip — ~gAB A BA"ip- 
By l(I)TJl ||B A || < || pa || 2^ < 0(A). Since also g& — > suffices to find a dense domain so that 

H-BaV'II = o(A _1 ). 

Now H.n can be thought of as symmetric functions in C 2 (M. 2N ). We take the subspace of 
functions in S(R 2Af ) which have a Fourier transform in QJ°(R 2Ar ) with support disjoint from the 
hypersurfaces where points coincide. If ip is in this space then a(pi)a(j>2)ip ~ 4'(PiiP2, ■ ■ ■ ) is a 
vector- valued function which has a Fourier transform in C^ D (M 2 x R 2 ) with support disjoint from 
the diagonal. Then 

u(p,q) = -jJ—a(^ + q)a(^-q)i; (92) 

has a Fourier transform u(X, x) which is an element of C^°(R 2 x (R 2 — {0})). Hence v = |x|~ 2 « 
is in the same space and if v(p, q) is the inverse Fourier transform then u = — A q v. 

Now we have for any n 

(B A ip){p) = J pA{q)u{p,q) dq 

= J pl{q)u(p,q) dq + 0(K- n ) (93) 

= j (-A q p* A (q))v(p,q) dq + 0{K- n ) 

Here we first replace the sharp cutoff pa by a smooth cutoff p* A and then integrate by parts. 
Since \A q p\{q)\ — 0(A~ 2 ) and since v(p,q) is rapidly decreasing in both variables we have 
\\B A ip\\ = C(A" 2 ) which suffices. 
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3. This follows by the Trotter-Kato theorem. 



Remarks. 

1. For N large our lower bound is H > — /i 2 e 167rAr2 . The coefficient 16ir can be improved but anyway 
the iV 2 behavior is probably not optimal. Indeed mean field calculations (Sj suggest that the 
actual lower bound may be e 0< - N \ The ground state is presumably a dense clump of particles: 
a "bosonic star" . 

2. For further studies of the spectrum on can consider the operator ^(E)^ 1 . We note that for 
E < if one scales all momenta by \J — E the operator becomes 



1 i r E 
-r- iog(— 

47T fi 



— log(—) + W (94) 



where 



W =(4vr)- 1 / dp X *(p) log (Ho + P 2 /i + l) x(p) 
1 



2(2tt; 



dpidp 2 dp 1 dp' 2 x*(pi +P2)x(p'i +P2) 



(95) 



a*(Pi)a* {p'2 ) 77 — ; ; ~ — —, — TTT a ( pi ) a ( p2 ) 

+*(pi -p'i)a*(p' 2 )-— ^ ,^ 1 a( P2 )) 

Ho + lji + uj 2 + u'. 2 + 1 / 

The issue is then to study properties of W. 
3.6 extensions 

We now allow an extra inter-particle potential v again assumed bounded. This means we add a 
potential 

V= - I a* (x)a* (y)v(x - y)a(x)a{y) dxdy (96) 



We have 

V X )\ < \\HooWNll\No - l) l / 2 M\\\Nl /2 {N, - 1)V2 X || (97) 

and thus V defines an operator on Hn satisfying ||V|| < A r2 |ju|| 00 /2. This is in configuration space 
and we actually consider the momentum space version V = r(^ r )V r r(J r_1 ) where T(T) is the induced 
Fourier transform on Fock space. This also satisfies ||V|| < A^ 2 ||t)|| 00 /2 which is the only fact we use. 
With a cutoff the full Hamiltonian is then 

H* = H + V' + H IA (98) 

Then Hf is self-adjoint on D(H Q ) n Hn- The same is true for 

H x = H + V (99) 

and we have H x > -7V 2 ||w||oo/2 

Proceeding as before we introduce resolvents 

R 1 {E) = (H 1 -EY 1 R*(E) = (H* - E)- 1 (100) 
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and for E < -N 2 ^]]^ 

$*(E) = gl 1 - B A R 1 (E)B* A (101) 

For such E we find as in lemma |S] that E is in the resolvent set of Hf on Ti, jy if and only if <&T (E) 
has a bounded inverse on TLn-2 hi which case 

R*{E) = + Ri{E)B\$t (E^BaR^E) (102) 



Theorem 4 



1. Let £7 < — e^- where 



e*=m^(N,N 2 \\v\\ 00 ,~^e 1 ^ N ^ c2 ^+^) (103) 

971. 

exists and is equal to 



and where C is the constant in lemma\^ For if; G TL^ the limit R&(E)ij) = liniA^oc R A (E)ip 



R*(E) = R^E) + R 1 (E)B*<I>*(Ey 1 BR 1 (E) (104) 

2. Rtt(E) is invertible. 

3. For E complex and not in [— e^,oo) the limit R#(E)ip = liniA-too B, A (E)ip exists. There is a 
self-adjoint operator iJ#(/j) with _ff#(/j) > — e% so that R&{E) = (iJ#(/i) — E)~ x . 

Proof. We follow the proof of theorem [3J We have 

R 1 (E) = R (E)-R 1 (E)V'Rq(E) (105) 

and hence 

= $ A (E) + B A R (E)V'R (E)B* A - B A R (E)V'Ri(E)V'Ro(E)B* A (106) 

For A = oo define $*(E) by replacing <f> A (E) by $(E) and B A by B. Since E < -^ e ^N 2 (c 2 \\v\u+i) 
we have for A sufficiently large or infinite instead of l|88|) 

$o,a(£) > 3N 2 (C 2 \\v\\oc + 1) (107) 

and it follows by the bound on $>i iA (E) that 

^a(E)>N 2 (C 2 \\v\\ 0O + 1) (108) 

For the other terms in I|1U6|) we note that E < —N implies \\B A R a (E)\\ < C by lemma |21 Also 
E < -N 2 |M|oo and the lower bound on H x imply that < (N 2 \\v\\ oc /2)- 1 . Using also 

||V|| < JV 2 |jz;|| 00 /2 

\\B A R Q (E)V'R (E)B* A \\ <C 2 ^1^/2 
\\B A R (E)V'R 1 (E)V'R (E)B* A \\ <C 2 N 2 \\v\\ OQ /2 
Combining these we find for A < oo 

so that <&a(£')~ 1 exists. Then for A < oo all E < —e* are in the resolvent set for R A (E) and (|102[> 
holds . 

As before $#(£") is essentially self-adjoint on VOHn On this domain $>#(E)ip = liniA^oo & A (E)%j). 
This follows from the convergence for <$> A (E) and the norm convergence of B a Rq(E). Using the uniform 
bounds on the inverses ^^(E) -1 converges strongly to ^(E)" 1 . 

Finally Rf(E) given by 11U2|> converges strongly to R A (E) given by 1)1041) . Here we use the norm 
convergence of B A R\(E) to BR\{E) which can be demonstrated using the adjoint of (|105fl . This 
completes the proof of the first part and the second and third parts follow as in theorem 



2Ar2 ( 109 ) 



<f>f(E)>N 2 (110) 
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A Trotter-Kato Theorem 



In the text we use the following version of the Trotter-Kato theorem . 

Theorem 5 Let E be a proper closed subset o/R and let H n be a sequence of self adjoint operators with 
resolvents R n (E) = (H n — E)^ 1 defined for all complex E ^ E. Suppose R n {E) converges strongly for 
some E ^ E and that the limit is invertible. Then there exists a self-adjoint operator H with resolvents 
R(E) = (H — E)^ 1 such that R n (E) converges strongly to R(E) for all complex E ^ E. 

A slightly different result is proved in 0. There E = R is allowed, but one needs convergence at 
two points with ±ImE > 0. This proof can be easily adapted to prove the quoted result. 
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